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ABSTRACT 


This program aims at developing mathematical models and computer codes 
based on these models, which allow prediction of the product distribution in 
chemical reactors for converting gaseous silicon compounds to condensed-phase 
silicon. The major interest is in collecting silicon as a liquid on the reac- 
tor walls and other collection surfaces. Two reactor systems are of major 
interest, the Westinghouse SiCl*/Na reactor in which Si(i) is collected on the 
flow tube reactor walls and the new AeroChem reactoi in which S 1 ( 2.) droplets 
formed by the SiCl*/Na reaction are collected by a Jet impingement method. 

During this quarter the following tasks have been accomplished: (i) par- 

ticle deposition routines have been added to the boundary layer code and 
(ii) Si droplet sizes in SiCl*/Na reactors at temperatures below the dew point 
of Si are being calculated. Initial results frou the latter show that droplet 
sizes in such reactors are predominantly in the range 0.2 to 1.0 uni. Such 
particles are too small to collect by conventional means (e.g., cyclones) but 
3 hould be easily collected by very high velocity impaction at reduced pressures 
such as that being examined by AeroChem. 
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1. IN TRODUCTION 


In thin report we consider the introduction of particles into the exist- 
ing model' for the vapor collection efficiency of silicon reactors of the 
West lnghouse tvpe. 4 The last report 1 concerned the partial collection onto 
the reactor walls of silicon vapor present in the reactor's hot core flow (see 
Fig. 1 for j schematic of the phvslcal situation!. While the deposition of 
vapor was attributed to the processes of convection, Fick diffusion, and 
turbulent (or eddy) transport to the reactor walls, the deposition of particles 
(silicon droplets' is expected to be controlled by additional mechanisms. In 
particular, one expects Brownian diffusion to govern the deposition of small 
particles ( 10"*-10“ * pm dlamt while Soret (or thermophoret ic' transport will 
become more important for larger particles (10“*-1.0 um diam', especially with- 
in the "viscous sublayer" region next to the reactor walls, where large tempera 
ture gradients exist. Away from the walls the particles will be transported 
entirely bv diffusive turbulent velocity fluctuations. This latter effect will 
always tend to drive the particles close to the walls. As was shown earlier,’* 
near the walls the former two molecular mechanisms (i.e.. Brownian and Soret 
effects' will predominate. 
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FIGURE 1 SCHEMATIC OF THE MODEL CONSIDERED 
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In order Co describe Che overall silicon collection efficiency (due to 
vapor and droplets) of a reactor, it is essential to model the deposition 
behavior of a distribution of particles’ varying in size from molecular dimen- 
sions to about a micron diameter. In the model considered here it is assumed 
that the particle size distribution function vPSDP for the silicon droplets 
entering the downstream section of the reactor for separation/collection is 
known.* The question, as before, is, "What fraction of the silicon vapor and 
droplets entering the reactor is collected on the walls?" The model described 
here attempts to answer this question, with due regard to the various deposi- 
tion mechanisms mentioned earlier and the two-dimensional developing boundary 
layer flow within the reactor. In fact, as a result of the slowly developing 
boundary layers within typical reactors most of the deposition occurs prior to 
the establichment of "fully-developed" conditions (i.e., with no streamwiso 
variations). Hence, the available simpler solutions which are valid for fully- 
developed, or partially developed (Graetz type)* flow cannot be used since such 
analyses only provide an asymptotic limit for the silicon mass transfer rate to 
the walls. Over the reactor's separation/collection length, the velocity, 
temperature, and concentration boundary layers are developing simultaneously. 
Hence, numerical solution of the coupled svstera of boundary layer conservation 
equations is unavoidable in determining the cumulative silicon collection effi- 
ciency, which depends upon the variation of the silicon mass flux (at the wall^ 
along the length of the reactor. 

This report outlines the general approach being adopted in the numerical 
modeling of the silicon separat ion/collection processes. Specific complications 
are isolated in order to show the nature of the problem and the methods used to 
overcome the difficulties encountered. The actual solutions for the collection 
efficiency will be presented in a subsequent quarterly report, which will 
discuss the model in greater detail. 

The results of calculations using the modified LAPP code are also discussed 
in this report. Particle size calculations using this code to describe a reactor 
similar to the Westinghouse reactor, but operating at temperatures below the Si 
dewpoint, show that the silicon droplets produced rapidly grow to a medium size 
*= 0.4 urn in diameter. 


* Actually, this PSDF will only be known after an analysis of the turbulent jet 
mixing, combustion, and nucleation processes in the upstream reactor section. 
This is, in fact, the main goal of the modified LAPP code being developed un- 
der this contract. 
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II. SILICON DROPLET DEPOSITION MODEL 

A, PHYSICAL AND MATHEMATICAL ASPECTS 


The present two-phase flow model for describing the transport of 
silicon particles from the hot ( « 3500 K) core-flow of the reactor to the 
cooled walls (* 1685 K) is based on the assumptions that: (i) the total parti- 

cle volume is considerably smaller th^n the volume of the surrounding gas phase 
sc that particle-particle Interaction is negligible; (ii> the concentration of 
the particles is small enough so that their presence does not alter the surround- 
ing gas flow field (i.e., velocity and temperature distributions); (ili) tie 
particles are small enough that they may be treated as "heavy-molecules" in 
local equilibrium with tiie gas. 

These assumptions are usually Justified in the downstream section of 
the silicon reactor helng considered. The dll. b eness of the silicon droplets 
can be expected to increase with Increasing distance from the droplet source, 
immediately following combustion (see Fig. 1). On the other hand, the nature 
of the processes leading to silicon droplet formation (gas-to-particle conver- 
sion via nucleation) ensures that the resulting particles will be of sub-micron 
sizes, in the absence of appreciable coagulation. 

Under these assumptions, one can write the following particle mass 
conservation equation suitable to the developing, boundary layer flow within 
the reactor (axisymmetric) : 
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Expressions for the Brownian diffusivity, D^, and the dimensionless Soret fac- 
tor, Oj, were presented earlier. * • k The turbulent diffusivity, Cp, is taken 
to be unaffected by the small particles being considered and hence given by the 
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expressions suggested by Lin et al. ? * Equation (1) ia valid pro- ided the 
particles do not grow appreciably (via vapor collection on their surfaces) dur- 
ing passage through the particle diffusion layer. Such an assumption may be 
justifiable in view of the typical thinness of these layers, as seen later. 
Moreover. Eq. (1) applies to particles of a uniform size, present in the 
external stream, depositing by convection, Brownian diffusion, and thermo- 
phoresis on the reactor walls. In order to describe the deposition of parti-.' 1 ', 
size distributions, therefore, one needs to repeatedly solve Eq. (1) with the 
appropriate boundary condition on concentration at the diffusion layer's edge, 
corresponding to each particle size class. That is, a PSDF expressing the 
variation of particle concentration vs. size in the external stream is dis- 
cretized into a finite number of independent size classes, with particles in a 
given class being characterized by uniquely defined transport coefficients, 
and D^. The velocity and temperature fields within the turbulent boundary 
layers on the reactor walls have to be supplied to the particle diffusion equa- 
tion. In the present case, for this purpose, the entire velocity and tempera- 
ture profile history along the reactor tube length is independently computed 
using a modification of the GENMIX code*’* for the finite-difference solution 
of coupled momentum and energy equations (Including variable fluid properties 
and the damping of turbulence near the reactor walls). 

Even in the absence of turbulence, the effects of thermophoresis and 
Brownian diffusion can immensely complicate the analysis of silicon droplet 
deposition behavior within the reactor. Physically, one expects thermophoretic 
effects to be important in regions of the reactor flow field where temperature 
gradients are large (i.e., near the walls). On the other hand, the larger 
particle size class, for which thermophoretic deposition is significant, corre- 
sponds to Brownian d if fusivities that are extremely small. There exists, there- 
fore, a thin sublayer region next to the walls where both Brownian diffusion 
and thermophoresis may be important for particle deposition. Outside this sub- 
layer, however, temperature gradients continue to remain sttong so that thermo- 
phoresis remains important, and diffusion by Brownian motion of particles 
becomes negligible. Thus one can consider the presence of an "inner layer," 


* This implies that the particles follow the turbulent eddy motions of the 
gas. Actually, particles near 1 um diam might break away from the eddies 
due to their even inertia. Such inertial effects are excluded in the pres- 
ent diffusional model. 
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next to the wills, where thermophoretlc transport of particles balances corvee 
tlon and Brownian diffusion. External to this is an "outer layer" wherein 
Cwiivcetlve transport balances thermophoretlc effects in conserving particle 
mass.* Although Eq . (1) is, in principle, general enough to describe both 
theae layers simultaneously, from the standpoint of computational efficiency 
it is useful to consider the abovementioned division into different layers 
(due to their disparate thicknesses). 

It has been established that, in analyzing the effects of thermo- 
phoresis, it is convenient to visualize the Introduction into the boundary 
layer of additional pseudo "sink" and "suction" behavior.* * The former will 
tend to deplete particle concentrations above the reactor walls while the 
latter acts towards increasing then. The net result of these opposing effects 
of thermophoresis is to appreciably increase the collection rate for the 
larger particles. 1 * 

B. CO M MENTS QN ..HIE— NUMER 1 CAL SOLUT I ON P ROCEDURE 

A few difficulties arise with the numerical solution of the particle 
transport problem, described in the previous section. They are briefly 
discussed in this section, together with strategies adopted to overcome them. 

1 . The Stiffness Problem 

When the characteristic time scales associated with convection. 
Brownian diffusion, and thermophoresis in Eq. (1) become widely disparate, one 
is likely to encounter "stiffness" during the numerical solution process. The 
problem is manifested by unstable numerical solutions if the integration does 
proceed or by a failure of the integrator to march forward, along the 
reactor length. The former corresponds to a build-up of truncation error as 
the solution proceeds (often seen is a wildly fluctuating solution) while the 


* It is to be noted that these layers art not related to the "viscous sublayer 
region next to the walls where turbulet.t dissipation of eddies occurs. In 
fact, both the "inner" and "outer" thermov'horetlc layers may be embedded 

within the "viscous sublayer" region. Moi ever, the division into "layers" 
adopted here is different from the asvmptoti' three-layer model implied by 
Refs. 10 and 11. 
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latter shows up as an inability of the integrator to meet the specified local 
error tolerance. 

To overcome this problem, a "stiffness solver," based on Gear's 
variable order, variable step-size, integration formulae is adopted. 1 * A 
method of lines approach is utilized, suitably discretizing the derivatives 
along the radial coordinate, lu order to solve Eq. (1) as an initial-value, 
ordinary differential equation. The "solver" employs an implicit scheme for 
greater stability and to allow marching in larger step sizes. Both these 
features are crucial to coping with the stiffness problem. 


The Scaling Problem 


Due to the disparate thicknesses of the "inner" and "outer" thermo- 
phoretic layers, discussed earlier, one runs into tabulation difficulties if 
the adopted radial grid {pacing approaches the thickness of the inner layer, 
or exceeds it. 


In order to deal with this problem the following transformation of 
the radial coordinate is defined: 



where r, r w , and r m are the radial locations corresponding to the point of 
interest, the reactor wall, and the edge of the particle diffusion layer. 

Note that the transformation restricts the domain of integration of Eq . (1) 
to the range 0^ r 1, regardless of particle size. It is, of course, neces- 
sary to specify the approximate thickness of the diffusion layer (i.e., 

5 ro “ r w “ r m ) • Since the thickness of the thermal boundary layer (i.e., ij'' , 
at a given x, is known via prior solution of the energy equation, one may 
readily estimate 6 n using the known Lewis number, for that particle size 
class, from the relation*: 

i si.(U )'/’ (5) 

m T g 


* A more cumbersome approach could have been to estimate 6 m by solving the 
boundary layer concentration integral equation 
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3 • Determ ination o f the Varia ble Edge Boundary Condi lion 

At th* above value of i n a concentration boundary condition needa to 
be specified. While the other boundary condition on Eq. (1) Is simply speci- 
fied at the wall as Y^ ■ 0, the one at the "inner layer" edge must be obtained 
via a degenerate form of Eq. (1). As discussed earlier, the degenerate equa- 
tion used here Ignores the diffusion term and furthe, assumes that the trans- 
verse convection of particles In the outer layer Is negligible compared to 
their axial convection. Thus one can write the "outer layer" equation as 


u 



K Y, 


( 6 ) 



( 7 ) 


and analytically solve it to obtain the following expresrion for the concentra- 
tion at .he outer edge of the "inner layer": 


Y 


i 



(83 


where the subscript 'o' refers to some specified Initial conditions (e.g., at 
tne inlet to the silicon separation/collect .on section of the reactor).* 
Equation (ft) emphasizes that the particle trajectories immediately outside 
the "inner layer" may be nearly parallel to the reactor walls. This is to 
be expected when the opposing thermophoretic and hydrodynamic forces acting 
on the particle nearly balance each other. 


* A more exact alternative approach to modeling the cuter layer will include 
the .lansverse convection term in Eq . (6) and, thus, require numerical solu- 
tion of differential equations for the concentration along particle 
trajectories. 


7 



III. SILICON DROPLET GROWTH RATES 


TN-207 


The acidified Lf~p code (now called CHEMPART) has been used to compute Si 
droplet formation and growth rates in en SiCl./Na reactor. The basic geometry 
is similar to that employed by Westlnghouse, * l.e., a 15 cm diam tube. Na 
vapor at 1400 K is input as a Jet of 7.5 cm diam with an outer flow of SiCl./ 
Hj/H/Ar at 3700 K. Input flows and conditions are summarized in Tuble 1. 

The mixing in the model flow is undoubtedly slower than in the actual Westing- 
house system since in the model SiCl. injection is tested as axial whereas in 
the actual reactor, SiCl. injection is transverse to the Na flow. A more seri- 
ous discrepancy may be that, to date, we have not input sufficient enthalpy to 
maintain the outer flow temperature at > 300 K. H and Hj immediately react 


TABLE I 

SiCl. REACTOR CALCULATIONS: INITIAL CONDITIONS 


Inner Flow 3 Outer Flow 3 


Velocity, cm s"‘ 

3550 

1642 

Temperature, K 

concentrations, mole fractions: 

1400 

3700 

SiCl. 

0.0 

0.381 

Na 

0.8 

0.0 

H s 

0.15 

0.429 

H 

0.0 

0.095 

Ar 

0.05 

0.095 


The inner flow is assumed uniform for the transverse 
distance r ■ 0 to 3.5 cm. The outer flow is uniform 
from r ■ 3.75 to 7.5 cm. Reactor radius ■ 7.5 cm, 
pressure - 1.0 atm, and the SiCl. to Na input flow 
ratio is stoichiometric. 


with SiCl. to form HC1 and SiCl 2 predominantly via the reaction scheme described 
in Table II. In doing ro the temperature rapidly drops to 2500 K. Equilibrium 
calculations show that, for the reactants and conditions of Table I, adiabatic 
reaction to equilibrium in the outer flow results in a temperature of 2350 K. 
Complete reaction with the Na in the inner flow leads to an adiabatic flame 
temperature of 2930 K. In future runs the enthalpy will be increased by increas- 
ing the H/Hj input ratio, and mixing of the jet and outer flow will be enhanced 
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Co more closely approximate Che accual reactor conditions. (A final flame tem- 
perature of *3200 K will be sought.) 

TABLE II 


REACTION MECHANISM 


1. 

Na ♦ 'iCl* 

* 

N«C1 ♦ SiCl , 

14. 

H ♦ SiCl* 

♦ 

HC1 ♦ 

SiCl, 

2. 

Na ♦ SiCl, 


NaCl + SiCl, 

15. 

H + SiCl, 


HC1 ♦ 

SiCl, 

3. 

Na + SiCl, 

•- 

NaCl ♦ SiCl 

lb. 

H ♦ SiCl, 

♦* 

HC1 + 

SiCl 

4. 

Na ♦ SiCl 

«♦ 

♦ 

NaCl ♦ Si 

17. 

H + SiCl 

«♦ 

•- 

HC1 + 

Si 

5. 

Si + Sx + M 


Si, ♦ M 

18. 

Si, ♦ Si, ♦ M 

-♦ 

Si, ♦ 

M 

6. 

Si ♦ Si a ♦ M 


Si, + M 

19. 

Si, ♦ Si, 

«♦ 

Si , , 


7. 

Cl ♦ Cl ♦ M 


Cl, + M 

20. 

Si,, + Si,, 


Si,* 


8. 

Na ♦ HC1 


NaCl + H 

21. 

Si, * + Si , * 


Si* , 


9. 

Na ♦ Cl a 

*- 

NaCl + Cl 

22. 

Si*, + Si*, 


Si,, 


10. 

H + HC1 


H, + Cl 

23. 

Si,, ■* Si,, 

-a 

Si , , , 


11. 

H + Cl, 


HC1 + Cl 

24. 

Si,,, + Si,,, 

♦- 

Si ,, * 


12. 

H + Cl + M 


HC1 + M 

25. 

Si,.* 

-* 

smallest particles 

13. 

H + H + ** 


H, + M 







At l ae lower (2500 K) temperatures now being calculated, Si particles form 
rrpidly and it is of interest to describe particle growth In such cases (even 
though, in the condensation mode employed by West inghouse , it is envisioned 
that they will operate above the Si dew point, i.e., at core flow temperatures 
well above 3200 K) , In Fig. 2 the temperature, and in Fig. 3, the Si atom mole 
fraction and Si droplet mass loading are plotted vs. flow time. The radial 
postion at which these plotted data were calculated is just beyond the edge of 
the inner flow, where reaction is initially the most intense. In Fig. 2 it is 
seen that in a very short distance, Reactions (7) and (12-17) of Table II result 
in a rapid temperature drop. The major species in the outer flow immediately 
become HC1 and SiCl,. As Na mixes into the region being described, the tempera- 
ture rises to 2750 K; incomplete reaction, and at this point, the absence of 
Si(i) prevents one from attaining the adiabatic flame temperature of 2930 K. 

As X/D increases beyond 1.0 wall losses (and radiation losses from particles to 
the 1800 K walls) begin to rapidly reduce the temperature. Figure 3 shows that 
at X/D ■ 0.8 the Si(g) concentration in the region bein t examined maximizes. 

The Si condensed pnase begins to form at this point and Si(g) is rapidly reduced 
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X/D 

llCURE 2 TEMPERATURE IN SiCl»/Na MIXINC REGION 


as the droplets grow via condensation. The calculation was carried out to the 

point at which about 85% of the available silicon, principally Si(g) and SiCl a , 

had been condensed. The droplets have rapidly grown, via condensation, to 

3ize = 0.25 urn. The mass and number densities (dm /d(ln r ) and dn /d(ln r ), 

P P P P 

respectively) at X / 1 ' - 2.8, the farthest downstream poin* - at which calculations 
were pursued, are shown in Fig. 4. From this point on, growth will proceed via 
agglomeration and will be slow. Roughly speaking, particle mass will double 
and number density will be halved for AX/D *■ 3000, i.e., particle growth is 
nearly over at X/D *» 5 and the size distributions of Fig. 4 are essential 1 '.’ fro- 
zen at this point. 


IV, CONCLUS ’ CNS 

A model presently being developed 1 * to describe the overall silicon collec- 
tion efficiency of reactors like the Westinghouse prototype has been outlined 
in this report. The model accounts for the deposition behavior of both silicon 
vapor and condensed liquid droplets via convection, thermophoresis, Brownian 
dif f usior., and turbulent (eddy) transport of particles. The model is general 
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FIGURE 4 Si DROPLET MASS AND NUMBER DENSITIES IN SiCU/Na 

MIXING REGION 

x/d - 2.8, r/d - 0.27 
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enough to handle the description of a particle size distribution function, pro- 
vided most of the particle growth has alreadv taken place prior to wall deposi- 
tion. Unlike existing op* roaches to this problem, ,0 » ‘ * the present model 
stresses the importance of several novel features that characterize a practical 
reactor. These includi a consideration of (1) variable fluid properties, 

(11) a description of the developing, turbulent boundary layers within the 
reactor, and (111) particle size distribution effects that invariable falsify 
deposition behavior,’ and hence the reactor's collection efficiency. 

In SiClfc/Na flow reactors In which input enthalpv is insufficient to 
achieve temperatures above the dew point of silicon, Si(4) droplet formation and 
growth calculations Indicate that droplet sizes will be subroicron. Most of the 
condensed Si in test cases now being examined resides in droplets of 0.2 to 
1.0 um diam. 


V, PLANS 

During the next quarter, boundary layer calculations will be made to more 
accurately determine Si collection efficiencies for the Westinghouse reactor, 
under a variety of operating conditions. The modified LAPP code (CHEMPART) 
will be used to predict Si droplet size distributions in the Westinghouse and 
AeroChem reactors so that this information can be used in the abovement ioned 
calculations. 


VI. NEW TECHNOLOG Y 


No reportable items of new technology have been identified. 
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